Interaction induced mergence of Dirac points in Non-Abelian optical lattices 
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We study the properties of an ultracold Fermi gas loaded in a square optical lattice and subjected 
to an external and classical non-Abelian gauge field. We calculate the energy spectrum of the 
system and show that the Dirac points in the energy spectrum will remain quite stable under on- 
site interaction of certain strength. Once the on-site interaction grows stronger than a critical value, 
the Dirac points will no longer be stable and merge into a single hybrid point. This mergence 
implies a quantum phase transition from a semimetallic phase to a band insulator. The on-site 
interaction between ultracold fermions could be conveniently controlled by Feshbach resonances 
in current experiments. We proposed that this remarkable interaction induced mergence of Dirac 
points may be observed in the ultracold fermi gas experiments. 
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One of the most interesting properties of graphene [l| , 
a single layer of carbon atoms packed in a honeycomb 
lattice, lie in the fact that the low energy excitations 
obey a linear dispersion relation [2J around the so-called 
Dirac points, and thus can be used as a testbed for the 
relativistic quantum electrodynamics. Consequently, it 
is now possible to observe many remarkable phenom- 
ena__in table-top experiments, such as Klein tunneling 
and the relativistic extension of Landau levels I5l- 
which usually only occur in high-energy physics [8|. 
This advantage of graphene has stimulated a great in- 
terest in the investigation of Dirac points Q in many 
other systems. In particular, ultracold atoms in optical 
lattices provide a versatile playground where the proper- 
ties of condensed matter systems can be simulated 0, [l(| 
in a highly controllable manner, such as the superfluid- 
Mott insulator transitions of Hubbard models |llMl3j . A 
quantum-optical analogue of graphene can be achieved 
by loading ultracold fermionic atoms such as 40 K or 6 Li 
lH, 13 1 into a hexagonal optical lattice [13]. The effects 
of Dirac points were discussed in the context of ulta- 
cold atoms in honeycomb lattice 14j and T3 (rhombic) 
lattices (l5| . Moreover, much more intriguing phenom- 
ena arise when these systems are subjected to artificial 
non-Abelian gauge fiel ds [ 16l 4l8l| , such as the non-Abelian 
Aharonov-Bohm effect ll 7| , non-Abelian atom optics [l9j . 
quasirelativistic effects [2fJ, or exotic topological phase 
transitions |21j . 

Here we would like to emphasize a fact that the non- 
Abelian artificial gauge field also provide an interesting 
setup where Dirac points emerge in a square optical lat- 
tice 12211 . which is originally limited to staggered fields 
[23I l24l ]. In this article, we consider a similar system 
in which two-component (two-color) ultracold fermionic 
atoms are trapped in a square optical lattice. And dra- 
matic difference between the two systems comes from the 
repulsive on-site interaction introduced into the model by 
us. Works 2l|, |22J mentioned above mainly dwelled on 
free Fermi gas on 2D optical lattice. We study the ef- 



fects of repulsive on-site interaction on the energy spec- 
trum of the 2D fermi gas loaded into a square optical 
lattice and subjected to a non-Abelian artificial gauge 
field. Implementing a self-consistent mean-field theory, 
we show that the Dirac points in the energy spectrum 
remain quite stable under repulsive on-site interaction 
of certain strength. When the on-site interaction grows 
stronger than a critical value, the Dirac points will no 
longer be stable and two Dirac points merge into a single 
one. This merging indicates a quantum phase transi- 
tion between a semimetallic phase and a band insulator 
[14 , l25l - l3l| . And one thing worthy to be mentioned here 
is that once the two Dirac points merge, the final dis- 
persion relation becomes quite exotic — it is linear in one 
direction but parabolic in the other orthogonal direction. 
Very recently, a well-designed experiment has been car- 
ried out by Leticia et al. [32J, in which the creating, mov- 
ing and merging of Dirac points has been realized with 
a Fermi gas loaded into a tunable honeycomb lattice. 
While the creating, moving and merging of Dirac points 
in the experiment 32| is generated by designing complex 
lattice geometries, the mergence of the Dirac points in 
our model is induced by strong repulsive on-site interac- 
tion. Since pairwise interactions could be conveniently 
controlled by means of Feshbach resonances [IH , we pro- 
pose that the interesting mergence of Dirac points in our 
model may be experimentlly observed and characterized 
in non-Abelian optical lattices. The non-Abelian opti- 
cal lattice could be prepared by generalizing the recent 
experiment 0, as proposed in 17l ljj]. 

We consider a two-component (two-color) Fermi gas 
trapped on a square optical lattice and subjected to an 
artificial non-Abelian gauge potential. Fermionic atoms 
in the system are interacting with repulsive on-site inter- 
action. It is well known that this pairwise interactions 
can be freely tuned by means of Feshbach resonances 33 1 
in nowadays ultacold atoms experiments. The Hamilto- 
nian of the system reads, 
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FIG. 1: (color online) (a) Two-component (two-color) ultra- 
cold fermionic atoms trapped in a square optical lattice and 
subjected to an non-Abelian artificial gauge field. The cir- 
cle filled with red(gray) color denotes internal atom state |t) 
U x {m) and U y (n + 1) are unitary operators induced 
by the external artificial gauge field, (b) The self-consisted 
mean-field order parameter A vs on-site interaction U. The 
system size is 24 x 24 and we set the convergence criterion as 
HT 4 . 
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where t is the hopping amplitude, c r ^(c\ a ) is the 
fermionic annihilation (creation) operator at site r of the 
square optical lattice, r can be regarded as pseu- 

dospin index, (rr'} denotes that the sum is over nearest 
neighbors and V is the strength of the on-site interaction 
between fermionic atoms. The coordinate of a fermion 
is given by r = (ma,na), where m,n are integers and a 
is the lattice constant of the square optical lattice. Here 
we set h = e = 1. The external gauge potential has 
the following form, A = ^-{—y,x) + a(B a a y , Bpa x ), in 
which Bo, B a , Bp are experimentally controllable param- 
eters and a X) y are the Pauli matrices. This intriguing 
artificial gauge field can be realized following the propos- 
als [121, |l8|, |34| , along the lines of the recent experiment 
[l6j |. After some algebra, the original Hamiltonian ([T]) 
becomes 



we constrain ourselves to the non-Abelian regime, where 
the Wilson loop \W\ < 2 and we set the Abelian flux 
$ = 0. 

In the noninteracting limit of Hamiltonian (JlJ, i.e. 
the case in which fermions hop freely between neigh- 
bor sites without any interaction, Hamiltonian of the 
system is a beautiful quadratic form and can be ana- 
lytically solved by Bogoliubov transformations. Energy 
spectrum of this case has been beautifully analyzed in 
literature 22], where the fermion gas becomes a col- 
lections of noninteracting quasiparticles and the spec- 
trum develops four independent Dirac points at Kjj £ 
{(0,0), (f ,0), (04), (f , I )} in the vicinity of marginally 
Abelian regime {§ a ,$p « tt/2). However, once the on- 
site interaction is taking into account in Hamiltonian ([!} , 
it's not a quadratic form any more and therefore can't 
be solved by Bogouliubov transformation directly. This 
is right the case we consider in this article. We study the 
repulsively interacting fermions on a square optical lat- 
tice subjected to a non-Abelian gauge field by the means 
of a self-consistent mean field theory. Our startpoint is 
Eq. ©. 

We mainly consider the repulsive interaction regime 
in this paper. As the on-site interaction grows stronger 
and stronger, fermionic atoms with different colors on the 
square lattice tend to repel each other and avoid staying 
on the same site. Once the interaction strength is over a 
critical point, the square optical lattice at half-filling will 
enter a phase in which each site of the lattice is occupied 
by single atom (see Fig. [TJ)). Therefore, we define A r = 

v{ c rt Cr -l-) as our order parameter. Under this mean- 
field approximation, the Hamiltonian ^ can be written 
as a quadratic form, 



H M F = - t E ^(ctAUrr^TT'Cr'T' + H.C.) 
<rr'> tt' 

- J^A^Crt + H.C.) + 
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(3) 



<rr / > r,r' 



(2) 

where U rr i is the matrix form of a nontrivial unitary oper- 
ator accompanying the hopping between nearest neigh- 
bor r and r'. If the hopping is along x axis, U rr > = 
U x (m) = e~ m97n e % ^ a<Ty . If the hopping is along y 

< a *. $ = B n a 2 is 



axis, U rr ' — U y (n) 
the Abelian magnetic flux, and 3> a ,/3 — B a ^a 2 is the 
non-Abelian flux. Fermions hopping around an elemen- 
tary square indicate a unitary transformation j22| U = 
U x (m)U y (n + l)U^.(m + l)U^(n) . The boundary between 
Abelian regime and non-Abelian regime is well defined by 
the gauge- invariant Wilson loop[22,[3ij W = trU. Here 



Through a canonical transformation, the above Hamilto- 
nian can be diagonalized by solving the following BdG 
equation [35[: 



E 



h rr ',f O rr ' 

0* r , 
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(4) 



where h Tr ',T = — t\Urr']rT, O rr ' = — A r 5 rr > — t[U TT ']ti 
and (m^/ , ) are the eigenvectors correspoinding to the 
eigenenergy E n , The self-consistent equation of the order 
parameter is 



Ar = F^u>"*tanh 



En 

2k B T 



(5) 
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We solve the set of BdG equations self-consitently via 
exact diagonalization method in real space. The system 
size of 24 x 24 is used in the calculation and the conver- 
gence criterion of A r is set to be in unit of nearest- 
neighbor hopping t. We find that the order parameter is 
uniform ( A r = A, where A is real) in the vicinity of the 
7r-flux regime. Our calculations (see Fig. [lj>) show that 
as the on-site interaction V increases from zero, the or- 
der parameter turn out to be non-zero at V c « 5.88i and 
the system undergoes a quantum phase transition from 
a semimetallic phase to a band insulator. 

By the above mentioned self-consistent mean-field 
theory, we transform the original Hamiltonian into 
a quadratic form Eq. (j3|). Implementing appropriate 
Fourier transformations, Eq. ([3]) can be easily diagnalized 
in momentum space. The corresponding energy spectra 
are shown in Fig. [2j As the mean-field order parameter A 
grows stronger, the two originally separate Dirac points 
(Fig. Ek) will first move closer(Fig. [2Jd), then merge into 
a single hybrid point at kh = (^,0). Around this hybrid 
point p — k — kh, the low-energy properties of the system 
are accurately described by a Dirac Hamiltonian 

H eff = J2 ^l H D^ P , H D = 2<r yPx - a xP 2 y (6) 
p 

where "fp = (c p f, c P i) T is the relativistic spinor. From 
this Hamiltonian we can see that the energy spectrum is 
linear in k x direction but quadratic in k y direction(Fig. 

The mergence of the two Dirac points signals a quan- 
tum phase transition from semimetallic phase to a band 
insulator [jy], [32| ■ If the order parameter A grows even 
stronger, a gap will be opened(Fig. which indicates 
a band insulator phase. The self-consistent mean-field 
theory calculations show that the order parameter A re- 
mains to be zero as long as the on-site interaction is 
smaller than U c — 5.88t. Once the strength of the on-site 
interaction is over U c , A will turn out to be a non-zero 
number which is not smaller than A c = 2. This can be 
easily seen in Fig. [T}d. Therefore, the spectrum shown in 
Fig. will not be observed in reality. We give out this 
spectrum in Fig. [2]just for comparison. 

To characterize the mergence of Dirac points more 
clearly, we calculate density of state(DOS) and cyclotron 
mass m c for different strength of on-site interaction (i.e. 
different values of A correspondingly). The results are 
shown in Fig. [3] and Fig. |4j For the density of states 
in Fig. |3l we find that dramatic difference appeared at 
the occasion when two Dirac points merge into a hybrid 
one. While the DOS displays two peaks for the case 
of A = 0(Fig. there are four peaks in the case of 

A = 2 (Fig. [3fc). Moreover, for A = 2 (Fig. Hi), DOS 
turns out to be non-zero at E = 0. These obvious differ- 
ences are supposed to be good indexes for the mergence 
of the Dirac points. In Fig. |4j we give out the cyclotron 
mass m c for different values of A, correspondingly, dif- 
ferent strengthes of on-site interaction. Comparing Fig. 
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FIG. 2: (Color online) Merging of Dirac points in the band 
structure of the system as A grows stronger($ a = &p = n/2). 
Top: Portrait of the energy spectrum in k x direction. Bottom: 
Variations process of the two Dirac points, (a) A = 0. There 
are two normal Dirac points, (b) A = 1.0. The two Dirac 
points moves closer, (c) A = 2.0. Two Dirac points merge 
into a single hybrid point, which signals the quantum phase 
transition, (d) A = 2.5. A gap is opened. 
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FIG. 3: Density of state (DOS) of the system vs. order pa- 
rameter A. (a) A = 0;(b) A = 1.0;(c) A = 2.0;(d) A = 2.5. 
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observed in ultracold fermi gas experiments nowadays. 
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and the NNSF of China (Grants No. 11075020 and No. 
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FIG. 4: (Color online) Cyclotron mass m c vs. A. (a) A = 
0;(b) A = 1.0; (c) A = 2.0;(d) A = 2.5. Inset: Comparison 
between (a) and (c) around the zero-energy point. 



and Fig. |4j;, we find that besides the difference in the 
numbers of peaks, the curve of m c changes from concave 
to convex around the zero-energy point as shown in the 
inset of Fig. Hfc. This change signals a quantum phase 
transition between the semimetalic phase and band insu- 
lator phase. 

From Fig. [TJd, we see that by increasing the strength 
of the on-site interaction, the mean- field order parameter 
can be tuned. And it's well-known that the on-site inter- 
action can be controlled conveniently by Feshbach reso- 
nances. That is, the mean-field order parameter could 
be modified by Feshbach resonaces indirectly. Therefore, 
this exotic mergence of the Dirac points induced by on- 
site interaction is supposed to be observed in ultracold 
fermi gas experiments. This is quite different from the 



works [311 . |32j , where they engineer Dirac points by tun- 
ing the nearest-neighbor tunneling [3lj | or by tuning the 
geometry of the optical lattices [32j. For the concrete re- 
alization of the experiment, one may use 40 K atoms in 
F = 9/2 or F = 7/2 hyperfine manifolds or 6 Li with 
F = 1/2. 

Dirac point plays a crucial role in many interesting 
phenomena in condensed-matter physics, for example, 
the massless electrons in graphene. In this article, we 
investigate a two-component (color) ultracold fermi gas 
which is loaded in a square optical lattice. The addition 
of non-Abelian artificial gauge filed gives rise of Dirac 
points in the energy spectrum of the system. We study 
the stability of Dirac points against the repulsive on-site 
interaction. Our calculations show that the Dirac points 
can be very stable under on-site interaction smaller than 
U c = 5.88i. At U Cl the Dirac points turn to be non-stable 
and merge into a hybrid point. The final hybrid point is 
linear in one direction but quadratic in the perpendicular 
direction. This mergence of the Dirac points denotes a 
quantum phase transition from semimetallic phase to a 
band insulator. This exotic phenomena is supposed to 
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